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ABSTRACT 

We present a study of the radial distribution of RR Lyrae variables, which present a range of 
photometric and pulsational properties, in the dwarf spheroidal galaxy Tucana. We find that the 
fainter RR Lyrae stars, having a shorter period, are more centrally concentrated than the more 
luminous, longer period RR Lyrae variables. Through comparison with the predictions of theoretical 
models of stellar evolution and stellar pulsation, we interpret the fainter RR Lyrae stars as a more 
metal-rich subsample. In addition, we show that they must be older than about 10 Gyr. Therefore, 
the metallicity gradient must have appeared very early on in the history of this galaxy. 
Subject headings: galaxies: dwarf — galaxies: individual (Tucana) — stars: horizontal-branch — 
stars: variables: other — Local Group 
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1. INTRODUCTION 

The oldest stellar populations in galaxies are of funda- 
mental importance because they can constrain the epoch 
of onset of star formation in the Universe. Unfortunately, 
the oldest main-sequence turn-offs (MSTO), which pro- 
vide the most straightforward characterization of the old 
population, can only be reached in a handful of Milky 
Way satellite galaxies using current ground-based observ- 
ing facilities. An alternative approach is the characteri- 
zation of the RR Lyrae star populations. The fact that 
the properties of the RR Lyrae stars reflect the properties 
of the old population is crucial because the horizontal- 
branch (HB) stars are about three magnitudes brighter 
than their MSTO counterpart. 

For example, it has long been suspected that the mean 
magnitude of the RR Lyrae variables, My(RR), is a 
function of [Fe/H], in the sense t hat more meta l-poor 
RR Lyrae stars are more luminous ()Sandagei fl958). The 
correlation between My(RR) and [Fe/H] was later con- 
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firmed by stellar evolution models, but with the addi- 
tional constraint that the morphology of the HB, i.e., 
the evolutionary status of t he stellar populati on, must 
be taken into account (e.g.. ICassisi et al.ll2004h . As the 
population gets older, the mass of the stars reaching the 
HB decreases, giving a bluer HB. In this case, the stars 
crossing the instability strip (IS) are therefore already 
evolving off the zero-age horizontal branch (ZAHB) and 
are more luminous than younger stars of the same metal- 
licity. 

In addition, observations of Galactic globular clus- 
ters suggested that the lower their metallicity, the 
longer the mean period of their RR Lyrae stars. This 
led to the various period - metallicity re l ations found in 
the literature (see, e .g., ICatelanl Il992t ISandaed 119931 : 
ISaraiedini et alj 120061 ) . or the mor e elaborate period- 
metallicity-am plitude relations (e.g.. lAlcock et all 120001 : 
Sandage 2004) , providing a convenient way of estimating 
the metallicity of a stellar system based on distance- and 
reddening-free observables (namely, mean period (P a b) 

and visual amplitude Avj. 

We note, though, that I Clement k Sheltonl (|1999t ) have 
argued that the period-amplitude relation for RRafr stars 
is not a function of metal abundance, but is instead re- 
lated to the evolutionary status of the RR Lyrae stars. In 
the metal-poor, Oosterhoff type II clusters, the (blue) HB 
stars would have started to evolve away from the ZAHB, 
crossing the IS at higher luminosities, therefore following 
a different PA relation fro m the ZAHB RR L yrae stars. 
However, theoretical dBono et alj 120071) and ob- 



serva ti onal (IWehlau et all 19991: iKaluznv et all 



| 2000t lLavden k Saraiedinil l2000t lOlech et all 120011: 
iBorissova et al. MM) evidence suggest that, while there 
is a clear correlation of Oosterhoff type with metallicity, 
there is an additional dependence of the period at 
a given amplitude upon metallicity, with the more 
metal-rich stars having shorter periods. 

An analysis of the variation of these observables as a 
function of the position within a galaxy can highlight the 
eventual radial trends across the studied galaxy. This, in 
turn, provides important clues to the formation mecha- 
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nisms and the star formation history of the host galaxy. 
Hence, by providing information about the properties 
of the underlying population, variable star research pro- 
cures a way to study the histories of these galaxies, in- 
dependent and complementary to the CMD analysis. 

In this Letter we report, for the first time, the detection 
of stellar population gradients in a dwarf galaxy based 
solely on the periods and luminosities of its RR Lyrae 
stars. This was made possible by the large number of RR 
Lyrae variables discovered in the isolated Local Group 
dSph Tucana, based on Hubble Space Telescope (HST) 
ACS data covering the whole body of the galaxy, together 
with the high quality of the light curves. We present sup- 
porting evidence from the different evolutionary phases 
present in the CMD. 

2. OBSERVATIONS AND DATA REDUCTION 

The present analysis is based on observations obtained 
with the Advanced Camera for Surveys onboard the HST 
as part of the major program LCID 11 aiming at deriving 
detailed star formation histories for a sample of Local 
Group isolated dwarf galaxies (C. Gallart et al. 2008, 
in preparation). As the goal of these observations was 
to reach the oldest main sequence turn-offs with good 
signal-to-noise on the final stacked images, 32 orbits were 
devoted to the observations of Tucana. These were col- 
lected over about five consecutive days between 2006 
April 25 and 30. As each orbit was devoted to one ex- 
posure in the F475W band and one in F814W, the ob- 
serving sequence consisted of alternating ~1000 seconds 
exposures in F475W and F814W for an optimal sam- 
pling of the light c urves. The DA OPHOT/ALLFRAME 
suite of programs ([Stetson! fl994) was used to obtain the 
instrumental photometry of the stars on the individual, 
non-drizzled images. A detailed description of the obser- 
vations, data reduction and calibration will be given in 
a forthcoming paper (M. Monelli et al. 2008, in prepa- 
ration). To sear ch for variable stars, th e Welch-Stetson 
variability index (jWelch &: Stetsonlll993[ ) was used. From 
the list of candidates, 358 stars were identified as bona 
fide RR Lyrae stars, for which the period was searched 
through Fourier analysis following the prescription of 
iHorne fc Baliunasl (|1986l ). We identified 216 RRa6, 82 
RRc, and 60 RRd based on their period and light-curve 
shape. Further details about the variable star content of 
Tucana and another isolated dSph, Cetus, will be pre- 
sented in a future paper (E. Bernard et al. 2008, in 
preparation) . 

3. THE RR LYRAE POPULATION AND ITS 

SPATIAL GRADIENTS 

Figure [1] ( top) shows the CMD of Tucana centered on 
the HB. The (F475W + F8UW)/2 ~ V filter combi- 
nation was chosen for the ordinate axis so that the HB 
appears approximately horizontal. The RR Lyrae stars 
are overplotted with larger dots using their intensity- 
weighted magnitudes. One can see that Tucana harbors a 
rather complex HB, which is well populated on both sides 
of the IS. The red side also presents a small gap in magni- 
tude, suggesting the combination of two HB of different 
luminosities. This is supported by the unusual width in 

11 Local Cosmology from Isolated Dwarfs, 

|http://www. iac.es/project/LCID/] 
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Fig. 1.— Zoom-in on the HB of Tucana. Top: The RR 
Lyrae stars are overplotted with larger dots. The solid black 
lines represent, from top to bottom, the Z=0.0001, 0.0003, 
and 0.001 ZAHB. The green and blue lines are evolution- 
ary tracks with Z=0.0003 and Z=0.001, respectively, and 
masses as labeled. Overplotted on them, the circles indicate 
intervals of 10 Myr. ZAHB lo ci and tracks are scale d-solar 
models from the BaSTI library l lPietrinferni et al.l l2004h , plot- 
ted assuming A„ =0.122 and A j =0.061 from NASA/IRSA 
(http://irsa.ipac.caltech.edu/applications/DUST/) and (m- 
M)o=24.84 (C. Gallart et al. 2008, in preparation). The 
vertical dashed lines roughly delimit the IS. Bottom: Different 
symbols have been used for each type and subsample of RR 
Lyrae variables, as labeled in the figure. For the sake of clarity, 
RRd are not represented. The horizontal line shows the mean 
intensity-weighted magnitude of all the RR Lyrae stars. 
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Fig. 2. — Period-amplitude diagram for the RR Lyrae stars in 
Tucana. Symbols are as in the bottom panel of Fig. [T] Linear fits 
to each sample of RRab are also represented. 



RR Lyrae stars in Tucana 



3 













All faint 






All bright 






Faint RRab 






Bright RRab 








Faint RRc+RRd -_ 








Bright RRc+RRd 




* 










\ 




















\ : 
















— "^^^^^^^^^ 















1 r 


f , r , . 1 , . 





50 100 150 200 

Semi-Major Axis [arcsec] 

Fig. 3. — Radial distribution of the different subsamples of RR 
Lyrae stars in Tucana, showing how the fainter RR Lyrae stars 
are systematically more concentrated in the center of the galaxy 
than the brighter variables. The error bars correspond to Poisson 
statistics. 

luminosity of the HB inside the IS (~0.3 mag at F475W- 
F814W=0.7, excluding the few bright and faint outliers). 
As the evolutionary tracks show, in Fig. [JJ (top), a stel- 
lar population at a given metallicity does not span more 
than ~0.15 mag in brightness within the instability strip, 
even taking into account the evolution off the ZAHB. 

We arbitrarily used the mean magnitude of the vari- 
ables at (F475W+F814W)/2. - 25.2, shown with the 
horizontal line in the bottom panel of Fig. [TJ to split the 
RRab, RRc and RRd variables into bright and faint sub- 
samples, each type of variable having approximately the 
same number of stars in each subsample. Different sym- 
bols have been used for each type and subsample. For 
the sake of clarity, RRd variables are not represented. 

The RR Lyrae stars are plotted on the period- 
amplitude (PA) diagram in Fig. [2] The symbols are 
the same as in the bottom panel of Fig. [TJ Linear fits 
to each sample of RRa6 are also represented. We can 
clearly identify two distinct RRa6 sequences, character- 
ized by different mean periods and dispersions around 
the fit. Based on their period distribution, we find a 
Kolmogorov-Smirnov (KS) discrepancy of Da's=0.4397, 
implying a probability Pks <10~ 9 % that the bright and 
faint RRab subsamples were drawn from the same parent 
population. As pre dicted by nonlinear pulsation models 
([Bono et al.lll997bl ). the mean period is a function of lu- 
minosity, in the sense that the more luminous variables 
tend to have a larger period. The period difference also 
show s up in the RRc su btype: the "bell-shape" distribu- 
tion ([Bono et al.lll997bh of the brighter RRc is shifted to- 
ward longer period with respect to the fainter RRc. The 
difference in mean period is ^0.02 days (Pks <0.34%). 
The mean primary period of double- mode RR Lyrae stars 
present a similar shift between the bright and faint sub- 
samples (P K s <2 x I CT 5 %). 

Interestingly, the faint and bright subsamples of RRa6 
variables have mean periods of 0.574 and 0.640 respec- 
tively, which is very reminiscent of the mean periods of 
the O osteroff types I and II (0.55 and 0.64, lOosterhofj 
Il939t ). This is also the case for the subsamples of RRc 
variables, having mean periods of 0.343 and 0.366, versus 
0.32 and 0.37 for the Ool and OoII globular clusters of 
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Fig. 4. — Color-magnitude diagrams of the inner and outer 
regions of Tucana, containing the same number of stars. The 
vertical lines delimit th e IS. Isochrones from the BaSTI library 
HPietrin fcrni et al. 2004), with ages and metallicities indicated in 
the figure, have been superimposed to highlight the gradients. Note 
the difference in the upper RGB, HB/red clump morphology and 
upper main-sequence. 

the Milky Way. On the other hand, the whole sets of 
RRab and RRc have mean periods of 0.604 and 0.353, 
respectively. These numbers would lead to the classifica- 
tion of Tucana as an Oost erhoff interme diate (Oo-Int), as 
is the case for most dSph (|Catclan 2004]) . In this particu- 
lar case, however, there is the possibility that the Oo-Int 
status could be attributed to the mixture of Ool and 
OoII populations, as hinted by the two roughly parallel 
sequences in the PA diagram of the RRab in Fig. O 

The presence of different populations in a galaxy, 
whether due to the details of its star formation history 
or to the accretion of an external stellar system, gen- 
erally leads to gradients in the observable properties of 
its stars. Figure [3] presents the radial profile for each 
subsample of RR Lyrae stars. The radii were chosen so 
that each concentric region contains the same number 
of variables. It shows that, for each type of RR Lyrae 
stars, the fainter variables are systematically more con- 
centrated near the center of the galaxy, while the brighter 
RR Lyrae stars are spatially extended. The difference in 
spatial distribution is very significant as the radii con- 
taining half of the bright and faint RR Lyrae stars are 
1.33' and 1.09', respectively, and a KS test gives a prob- 
ability of <0.001 for the bright and faint RR Lyrae stars 
to have the same radial distribution. The combination of 
different intrinsic properties of the individual stars with 
the different spatial distribution supports the hypothesis 
that they represent separate populations. 
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Fig. 5. — Mean mass as a function of age and metallicity for 
artificial stars within the IS of Tucana, from IAC-STAR simula tions 
using the BaSTI library with r;=0.4 (Pictrinfcrni ct al. 2003) and 
no dispersion in mass loss at a given age. 



4. DISCUSSION 

The difference in the mean period of the RRa6 stars 
of each sample could be attributed to both a difference 
in metallicity (from the period-metallicity correlation) or 
a difference in the evolutionary status of the individual 
stars, the stars in the redward evolution off the ZAHB 
having a longer period (jBono et al. 1997b], their Fig. 16- 
17). However, the top panel of Fig. [1] shows that evolu- 
tion off the ZAHB alone cannot account for the range of 
luminosity spanned by the RR Lyrae stars. Hence, it is 
necessary to invoke a range of metallicity to reproduce 
the distribution of stars within the IS. The hypothesis 
of a bimodal metallicity distribution is strengthened by 
the double "bell-sh ape" of the RRc in the PA diagram 
(|Bono et al.lll997ah and the presence of two RGB bumps 
separated by ~0.2 mag in F814W (C. Gallart et al. 2008, 
in preparation) . We use the luminosity of the theoretical 
ZAHB in Fig. Q] to estimate this range. A firm upper 
limit for the metallicity of the old population is given by 
the ZAHB at Z=0.001. While the Z=0.0001 ZAHB seems 
to be slightly too luminous, the Z=0.0003 tracks cannot 
produce the brightest RR Lyrae stars in our observations. 
We suggest Z~0.0002 as an approximate lower limit to 
the metallicity of Tucana. 

Although the range of metallicity is relatively wide for 
a dwarf spheroidal galaxy, it is supported by the color of 
the RGB of the inner and outer regions of Tucana pre- 
sented in Fig. |U The isochrones show that the higher 
metallicity population visible in the central region, with 
Z~0.001, is absent in the outer region, while a metallic- 
ity slightly lower than Z=0.0003 is needed to explain the 
blue side of the sub-HB RGB. Note, also, the difference 
in the morphology of the HB. The faint RR Lyrae vari- 
ables and the red HB stars show up only in the center 
of the galaxy, while the brighter stars of the IS and the 
blue HB stars are present everywhere. Consequently, the 
HB ratio (HBR 12 ) increases by about 0.5 over the radius 
sampled with our data, similar to what was observed 
in the Sculptor dSph dMaiewski et al.l I1999L see also 
iHurlev-Keller et all 1 19991) for which spectroscopy indi- 



cates the presence of a metallicit y gradient (iTolstoy et all 
l2004h and a metallicity spread (jClementini et al.ll2005h 
within the old (>10 Gyr) HB population. 

Our metall ic ity r ange supports the results of 
lHarbeck et all (|2001l ). who found a bimodal [Fe/H] 
distribution in Tucana from the RGB stars over 
their observed HST WFPC2 field. From isochrones 
in the HB morphology-metallicity plane, they ex- 
plain the gradient in HB morphology as a 'pure 
metallicity effect'. Howe ver, artificial CMD s pro- 
duced by IAC-STAR flA paricio & Gallart] l200i 
|http://iac-star.iac.es/iac-star/ ) show that to obtain 
a red HB component such as the one observed in the 
center of Tucana, stars having a total mass in the 
range ^0.75-0.80 M Q are needed (see, also, the top 
panel of Fig. [T]). Assuming the common value for the 
mass loss parameter 77= 0.40 with the BaSTI library 
(|Pietrinferni et ail I2004T ) . these masses correspond to 
stars younger than about 9 Gyr old, quite independently 
of the metallicity between Z=0.0001 and 0.001. The 
same models give an age of ~13 Gyr for the bluest HB 
stars at Z=0.0001 which are present over the whole 
galaxy. 

The morphology of the MSTO in Fig. Q] reflects this 
range of age: one can see that the CMD of the central 
region harbors a younger population than the outer re- 
gion. Fitting isochrones to the subgiant branch gives 
a population age ranging from about 11 to 14 Gyr old 
in the outer part of the galaxy, while it appears to be as 
young as ~8 Gyr in the center. Note that crowding tends 
to broaden the features in a CMD and may be partly de- 
creasing the estimated age of the youngest MSTO stars 
in the center of the galaxy. Nonetheless, a preliminary 
SFH of the inner and outer regions of Tucana — which 
takes into account the effects of crowding — do show that 
there is a slightly younger population in the center. 

However, the youngest stars within the IS are slightly 
older than the red HB stars. Figure [S] illustrates the re- 
sults of our artificial HB computations: it represents the 
mean mass as a function of age and metallicity for arti- 
ficial stars within the instability strip of Tucana, and in- 
dicates that the faintest, highest metallicity (Z>0.0006) 
variables of Tucana must be >10 Gyr old. Therefore, 
under the reasonable assumption that chemical enrich- 
ment follows age in star forming galaxies, the presence 
of gradients in the RR Lyrae populations shows that these 
metallicity gradients appeared very early on in the history 
of this galaxy. 

Facility: HST (ACS) 
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12 HBR = (B-R)/(B+V+R), where B, V, and R are the numb ers 
of stars to the blue, within, and to the red of the IS (Lee 1990). 
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